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ABSTRACT 

We present the 1.9-4.2 /im spectra of the five bright (L < 11.2) young stars 
associated with silhouette disks with moderate to high inclination angle of 39-80° 
in the M42 and M43 regions. The water ice absorption is seen toward dl2 1-1925 
and d216-0939, while the spectra of dl82-316, dl83-405, and d218-354 show no 
water ice feature around 3.1 /im within the detection limits. By comparing the 
water ice features toward nearby stars, we find that the water ice absorption 
toward dl21-1925 and d216-0939 most likely originates from the foreground ma- 
terial and the surrounding disk, respectively. The angle of the disk inclination 
is found to be mainly responsible for the difference of the optical depth of the 
water ice among the five young stars. Our results suggest that there is a critical 
inclination angle between 65° and 75° for the circumstellar disk where the water 
ice absorption becomes strong. The average density at the disk surface of d216- 
0939 was found to be 6.38 x 10~^^ g cm~^. The water ice absorption band in the 
d216-0939 disk is remarkable in that the maximum optical depth of the water 
ice band is at a longer wavelength than detected before. It indicates that the 
primary carrier of the feature is purely crystallized water ice at the surface of the 
d216-0939 disk with characteristic size of ~ 0.8 /im, which suggests grain growth. 
This is the first direct detection of purely crystallized water ice in a silhouette 
disk. 
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1. INTRODUCTION 

Planetary systems are thought to evolve from circumstellar disks around young stellar 
objects. In this context, dynamical and chemical evolution of circumstellar disks has been 
extensively studied with various observational and theoretical methods to reveal how plan- 
ets are formed. Recent advanced astronomical facilities of both ground-based 8- 10m class 
telescopes and space missions such as Hubble Space Telescope (HST), ISO, AKARI, and 
SPITZER have revealed morphological and spectroscopic features of circumstellar disks in 
great detail with high sensitivity and high spatial resolution. 

The most striking discovery was the direct detection of circumstellar disks in Orion 
Nebula Cluster imaged by HST in its early operation phase (O'Dell et al. 1993). They 
exhibited clear morphology of the circumstellar disks as a silhouette in front of the very bright 
background nebula illuminated by O, B stars (Trapezium stars) in the cluster. This was a 
huge milestone for studies of the protoplanetary disks, as it showed that the circumstellar 
disks could be directly studied with high spatial resolution (~ O'.'l). The silhouette disks 
have definite advantages, because various basic properties of the circumstellar disk such as its 
shape, size, inclination angle, and composition are unambiguously determined. It provides 
us with excellent opportunities to investigate the circumstellar disk as a progenitor of the 
planetary system. 

Water is a key element at the top of the hierarchy of building blocks to complex, pre- 
biotic molecules. In our Solar System, water is thought to be brought onto the surface of a 
terrestrial planet such as the Earth in the form of water ice in its sohd phase (van Boekel 
2007). Therefore, the search for the water ice in circumstellar disks is of extreme interest to 
investigate the possible emergence of the hfe elsewhere. The water ice feature at the short- 
est wavelength appears at 1.5 /xm, 2.0 //m, and 3.1/xm in the near infrared region, which is 
accessible from the ground. In particular, among those features, the water ice absorption at 
3.1 /xm is the most important in terms of its detectability. 

The water ice absorption at 3.1 //m has been detected toward various astronomical 
sources including young stellar objects since its first detection by Gillett & Forrest (1973). 
While in most of the cases the detected water ice is attributed to the foreground cloud and/or 
their envelope around the sources, recently the water ice absorption in the circumstellar disks 
has been detected successfully (Pontoppidan et al. 2005; Terada et al. 2007; Schegerer & Wolf 
2010; Aikawa et al. 2012). However, those detections are based on the scattered light by the 
circumstellar particles around the central star, and the geometry of the water ice absorption 
is not clear. Water ice absorption seen directly through the silhouette disks will be easier 
to interpret than those detections. Another advantage comes from the fact that the M42 
and M43 regions suffer from less extinction (Scandariato et al. 2011), which means that 
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the foreground contamination is minimized for the measurement of the water ice in the 
circumstellar disk. Taking into account the fact that majority of stars are born in high mass 
star forming regions, it is important to investigate the closest high mass star forming regions 
of M42 and M43. The silhouette disks were initially discovered in the Orion Nebular Cluster; 
M42 (McCaughrean & O'Dell 1996; Bally et al. 2000), and recently the sample is largely 
expanded with several additional silhouette disks found also in M43 (Smith et al. 2005). In 
this paper, the distance to M42 is adopted as 414 pc (Menten et al. 2007), while 436 pc is 
used for M43 (O'Dell & Henney 2008). 

Wc present the spectra of five sources with silhouette disks in the M42 and M43 regions, 
namely dl21-1925, dl82-332, dl83-405, d218-354 (Bally et al. 2000), and d216-0939 (Smith 
et al. 2005). In §2, the observational procedures and applied data analysis are described. In 
§3, the continuum of each target is determined and the optical depth of the water ice feature 
is derived. The origin of the detected water ice, the possible cause of the different optical 
depth of the water ice, and the interesting nature of the detected water ice in d216-0939 disk 
are discussed in §4. The conclusion is described in §5. 

2. OBSERVATION & DATA REDUCTION 

The combination of near infrared imaging (1.4-4.2 /xm) and spectroscopy (1.9-4.2 /xm) 
were conducted as a series of the observations from 2000 December to 2007 September. All 
the observations were performed under seeing limited conditions using the Infrared Camera 
and Spectrograph (IRCS; Tokunaga et al. 1998; Kobayashi et al. 2000; Terada et al. 2004) 
at the Subaru 8.2m Telescope atop Maunakea, Hawaii. IRCS was initially mounted at the 
Cassegrain focus of the telescope in 2000, and moved onto the Nasmyth focus on 2005 
December for the upgraded adaptive optics system (A0188; Hayano et al. 2010). The pixel 
scales became slightly finer by a factor of 0.89 after moving to the Nasmyth focus. The 
observation log is shown in Table 1 with a seeing size on those nights. The seeing values 
were good in most of the observing runs. But the seeing was poor in some cases, especially 
for d218-354. Sky conditions throughout all of the observing nights were photometric. 

Imaging data at L' were obtained for all the objects, and in addition imaging at H 
and K were carried out for d2 16-0939. For all the imaging, a box-shaped five- image dither 
pattern was used with a separation of 5'.'0 and sky frames were obtained 300"west away from 
the objects. Spectroscopy for dl21-1925, dl82-332, dl83-405, and d218-354 was conducted 

at the Cassegrain focus with a slit width of 0'.'6, corresponding to a spectral resolution of 
A/AA ~ 360 and 220 at K and L, respectively. Only for d216-0939, the slit width of 0'.'45 
was used at the Nasmyth focus to obtain the spectral resolution of 430 and 260 at K and L. 
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Table 1. Observing Log 





Observing 




Exposure 


Seeing 


Standard 


Object 


Date 


Mode 


Time 


Size 


Star 




(UT) 




(sec) 


(") 




dl21-1925 


2001 Mar 02 


L' Imaging 


90 


0.58 


HD 36512 




2001 Mar 10 


K Spectroscopy 


480 


0.96 


HR 1931 




2000 Dec 04 


L Spectroscopy 


2520 


0.45 


HR1931 


dl82-332 


2001 Mar 02 


L' Imaging 


90 


0.66 


HD 36512 




2004 Feb 04 


K Spectroscopy 


480 


0.19 


HR 1931 




2001 Jan 11 


L Spectroscopy 


1250 


0.31 


HR 1931 


dl83-405 


2001 Mar 11 


h' Imaging 


765 


0.91 


HD 36512 




2004 Feb 04 


K Spectroscopy 


480 


0.32 


HR 1931 




2001 Jan 11 


L Spectroscopy 


1200 


0.66 


HR 1931 


d2 16-0939 


2006 Jan 15 


H Imaging 


90 


0.49 


FS 13 




2006 Jan 15 


K Imaging 


90 


0.52 


FS 13 




2006 Jan 15 


L' Imaging 


180 


0.42 


HD 40335 




2006 Jan 15 


K Spectroscopy 


480 


0.44 


HR 2328 




2006 Jan 15 


L Spectroscopy 


2160 


0.42 


HR 2328 


d2 18-354 


2001 Mar 11 


L' Imaging 


405 


1.29 


HD 36512 




2001 Mar 10 


K Spectroscopy 


480 


0.81 


HR 1931 




2001 Jan 10 


L Spectroscopy 


1400 


0.90 


HR 1931 



Comparison Observation for Vicinity Stars 



JW 370 2007 Sep. 06 V Imaging 135 0.27 HD 225023 

2007 Sep. 06 L Spectroscopy 1260 0.27 HR 1724 
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The position angle of the sht was set to the perpendicular direction of the silhouette disk for 
dl21-1925, dl82-332, dl83-405, and d218-354. In the case of d216-0939, the position angle of 
the slit was chosen to obtain the nearby star JW 671 (Jones & Walker 1988) simultaneously. 
Standard A-BB-A sequence was applied to subtract the sky emission with an offset width of 
7^0. 

Regarding flat-ficlding, imaging at V used a sky flat. A flat field image at H and K 
and also flat field spectrum at K and L were taken with an integration sphere illuminated 
by a halogen lamp. 

In addition two stars, JW 370 (Jones & Walker 1988) and JW 671, which are located 
at the vicinity around dl21-1925 and d216-0939 respectively are observed for comparison. 
The observing parameters for JW 370 are included in Table 1. The spectrum of JW 671 was 
obtained simultaneously with d216-0939 on the same slit. 

Standard data reduction processing was applied using the IRAF software packages for 
flat fielding, sky subtraction, telluric correction, and wavelength calibration. The standard 
stars for telluric correction are listed in Table 1. Applied photometric values at V for HD 
36512, HD40335, and HD 225023 arc 5.54, 6.44, and 6.979 mag, respectively (Leggett et al. 
2003). The magnitudes at H and K for FS 13 are assumed to be 10.176 and 10.126 (Casali & 
Hawarden 1992). For spectroscopic reference, the spectral type of HR 1931, HR 2328, and HR 
1724 are 09.5V, AOVn, and AOV, whose spectra are assumed to have a blackbody radiation 
with a temperature of 34,600K, 9,480K, and 9,480K, respectively. The wavelength cahbration 
was performed using the telluric absorption lines in the spectrum. Since spectral region 
around hydrogen recombination lines (Br5; 1.94508 //m, Br7; 2.16611 /xm, Pfe: 3.03919 /im, 
Pf5; 3.29698 /xm, Pf7; 3.74054 /xm. Bra; 4.05224 yum) shows insufficient sky subtraction due to 
high spatial variation of the strong nebular emission from the H II region, those wavelengths 
arc omitted from the spectra except for d2 16-0939 which has rather weak nebular emission. 
The hydrogen absorption of spectroscopic standard star HR 2328 is appropriately removed 
using the method developed by Vacca et al. (2003). The differences of the airmass between 
the target and the standard for all the spectroscopy were kept to be less than 0.04, while 
imaging data were taken with an airmass difference of < 0.13. 



3. RESULTS 

For imaging data, the aperture photometry was performed with a radius of 1'.'5 except 
for the case of d218-354 at V with a radius of 2'.'5. The sky region for the determination of 
zero level was selected carefully for each target because the observed region is nebulous and 
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confused by many sources. It is essential to have a good image size under < I'/O. Thanks 
to the good seeing conditions almost throughout the observing runs, the confusion from the 
nearby stars did not cause serious problems. The only exception is for the case of d218- 
354. The seeing condition was poor for the observation of this target, which caused worse 
quality of the data especially at L' . Since the area of d218-354 is less confused by the nearby 
stars than the case of the others, there was no problem even in such bad seeing conditions. 
The photometric results are summarized in Table 2 together with some of the values from 
the literature. Our infrared broad-band images do not show the silhouette disks, which are 
seen in the Ha narrow-band images taken by HST. Therefore, the silhouette disks have 
a neghgible effect on our photometric results of the continuum of the central sources. The 
V images of dl21-1925, dl82-332, dl83-405, and d218-354 show a point source profile, and 
the fluxes for those objects originate from the central source. Also in the case of d216-0939, 
which exhibits distinct scattered light in the HST images, our infrared images at H , K, and 
L' do not clearly show the morphology of the circumstellar disk. The scattered light of the 
d216-0939 images is discussed in more detail in §4. 

All the extracted 1.9-4.2 /im spectra for the targets associated with the silhouette disk 
are shown with black lines of Figure 1. These spectra in the atmospheric windows of K-band 
(1.9-2.5 /xm) and L-band (2.8-4.2 /xm) were adequately offset to match to the photometric 
data determined in the above. Accurate photometric values at K and L' are very critical 
to precisely extract the possible absorption feature around 3.1 iim region. For dl21-1925, 
dl82-332, dl83-405, and d218-354, the photometry at K taken at the different epoch is 
used, and the possible variability of each target should be carefully taken into account. 
Regarding K magnitudes for dl82-332, dl83-405, and d218-354, McCaughrean & O'Dell 
(1996) and Muench et al. (2002) shows consistent values, which suggests those are not 
variable at K. As for dl21-1925, there is no photometric data available at K other than 
the one by McCaughrean & O'Dell (1996) to investigate the possible variability. However, 
our spectroscopic continuum level shows sufficient continuity and it should be acceptable 
to assume that the K magnitude measured by McCaughrean & O'Dell (1996) is useful for 
the following estimate of the absorption feature. The derived continuum is confirmed to 
be in good agreement with the available photometry at H except for d2 16-0939, which is 
brighter by 0.3 mag most probably due to contribution of the scattered light around the 
central source. 

In Figure 1, the broad absorption feature around 3.1 /im is definitely seen in the spectrum 
of d216-0939. For the other sources, the detailed determination of their continuum emission 
is obviously required. To determine the continuum level, the spectral template in the IRTF 
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Fig. 1. — 1.9-4.2 /iin spectra of five silhouette disks (black lines). The model spectra for 
fitting for the spectra of the targets are shown with gray lines. d2 16-0939 clearly demon- 
strates the broad absorption band around 3.1 fim, which is attributed to the water ice. The 
same feature can be marginally seen in the spectrum of dl21-1925. The right panel shows 
the HST images taken with a narrow band filter of F656N at Ha line. Images dl2 1-1925, 
dl82-332, dl83-405, and d218-354 are from Bally et al. (2000), and the d216-0939 image is 
from Smith et al. (2005). All the figures are displayed in the same scale of 2'.'73 x 2'.'73, and 
the scale bar of 0'.'5 is shown inside each image. 
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Spectral Library (Rayner et al. 2009) which is available on-line ^ is adopted. For dl21- 
1925, dl83-405, and d218-354, the spectral type and the visual extinction (Ay) which are 
summarized by Getman et al. (2005) from the original data (Hillenbrand 1997; Luhman 
et al. 2000) are applied as an initial assumption. The spectral type and the Ay for d216- 
0939 are estimated by Hillenbrand (1997). All this information is summarized in Table 3. 
Meanwhile, no reliable spectral type is found for dl82-332, and only Ay by McCaughrean 
& O'Dell (1996) is presented in the table for this object. To obtain the good continuum 
fitting, it is well known that some contribution of the infrared excess/ veiling effect should be 
taken into account. Here, the single temperature blackbody radiation is added to reproduce 
the measured spectral continuum. This value represents the characteristic temperature of 
the material attributed to the longer wavelength continuum. In fitting the continuum, the 
parameters for the already known spectral type and Ay were tried first. For dl21-1925, 
d218-354, and JW 671, the spectra can not be fit using the already known parameters, and 
those were adjusted to obtain a good fit. The best fit parameters for the spectral model are 
shown in Table 3, and the derived optical depths of the water ice are presented in Table 4. 
Regarding the spectral type and Ay, the parameters that differ from previous works are 
represented with the bold font in Table 3. The detailed spectral classification is not covered 
by this paper, especially because the classification in the infrared is usually less accurate 
than in the optical. Since the discrepancy is small enough, it will not affect the following 
overall discussion. 

For comparison in the discussion section below, the spectra of two stars: JW 370 and 
JW 671 are shown in Figure 2 with an overplot of the estimated model spectra. Exactly 
the same procedure is applied for determination of the possible absorption around 3.1 /im 
region. Since no spectral information in K for JW 370 is obtained, the photometry data at 
K is used for JW 370. The spectral type and the Ay for those stars derived by Hillenbrand 
(1997) can be found in Table 3. 



4. DISCUSSION 

4.1. Origin of Detected Water Ice toward dl21-1925 and d216-0939 

dl2 1-1925 is located in one of the most embedded region [Ay > 6) in the Orion Nebula 
cluster (Scandariato et al. 2011), and contributions of the water ice material in the foreground 
to the detected water ice should be carefully evaluated. In fact, the shallow, but definite 



^http : //irtfweb.ifa.hawaii.edu/^spex/IRTF.SpectralJLibrary/ 
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Fig. 2. — Spectra of JW 370 and JW 671. These are shown in the same range of 1.9-4.2 /im 
for comparison. JW 370 and JW 671 are the vicinity stars of dl21-1925 and d216-0939, 
respectively. The black filled square is the K band flux measured from acquisition images. 
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water ice absorption at 3.1 //m was seen in the spectrum of JW 370 which is located only 
3'.'39 west (1403 AU in projected distance) away from dl21-1925 in Figure 2. For comparison 
between those two detections of the water ice, Figure 3 shows the optical depths of the 
detected water ice at 3.1 fira for both dl21-1925 and JW 370. The profile and the depth of 
those two features are found to be matched well within the signal-to-noise ratio, and it is 
most likely that the detected water ice in both dl2 1-1925 and JW 370 could be attributed 
to the common origin of the foreground material. 

In contrast, the Scandariato et al. (2011) data shows less extinction (Ay ~ 3) in the 
region where d216-0939 resides. The nearby star JW 671 is 4'.'38 south-west (1910 AU in 
projected distance) away from d2 16-0939 and is used for evaluation of possible extension of 
the water ice absorbtion around the d216-0939 disk. Figure 2 shows that JW 671 has no 
water ice absorption, and the difference of the optical depth of the water ice in d216-0939 
and JW 671 is presented in Figure 3. JW 671 is known to be a member of the Orion Nebula 
cluster (Furesz et al. 2008), and thus it is not a foreground star. Therefore, the water ice 
detected toward d216-0939 is concluded to be locahzed at least within a scale of 2000 AU 
from the central star, which is comparable to a typical size of the disk-envelope around young 
stars. 



4.2. Possible Explanation for Water Ice Detection in d216-0939 Disk 

Silhouette disks in M42 and M43 represent a wide variety of characteristics in its mass, 
size, and inclination angle (McCaughrean & O'Dell 1996; Smith et al. 2005; Bally et al. 
2000; Ricci et al. 2008; Mann & Williams 2010). In addition, the central stars associated 
with the silhouette disks also have range of spectral type, mass, and Ay (Hillenbrand 1997; 
Luhman et al. 2000) as well as various distance from the O, B stars (Trapezium cluster; 
mostly 9'^ Ori C for M42 and NU Ori for M43) exciting the H II regions (Vicente & Alves 
2005). Parameters of the disk mass, size, inchnation angle, and distance from the exciting 
stars are summarized in Table 5. In the following, the relationship of the water ice detection 
with each parameter is discussed. 

4.2.1. Disk Mass 

Mann & WiUiams (2009) derived the disk mass of (4.50 ± 0.06) x 10"^ Mq for d216-0939 
from the submm measurement at 880 /im. This is the second most massive circumstellar disk 
among their large sample of protoplanetary disks in M42 and M43, and the most massive in 
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Fig. 3. — Comparison of the optical depths of the water ice toward the objects (dl21- 
1925 and d216-0939) and the vicinity stars (JW 370 and JW 671, respectively). Since the 
depth and the profile are very similar between dl21-1925 and JW 370, a common origin 
from foreground material is suggested for both water ice. On the other hand, the features 
are definitely different for the case of d216-0939 and JW 671. Localized water ice in the 
circumstellar disk of d216-0939 is responsible for the difference. 
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our five targets. In addition, the d2 16-0939 disk has the second largest size among all the 
silhouette disks in M42 and M43 (Vicente & Alves 2005; Smith et al. 2005). It might be 
straightforward to expect a larger column density of water ice in the bigger disk. However, 
taking into account the null detection of the water ice in the d218-354 disk which is the 
second most massive and the second largest in our sample, this is not a unique reason for 
finding water ice only in the d2 16-0939 disk. 

4-2.2. Distance from O, B stars 

The distance to exciting sources for H II regions is a key parameter for determining the 
water ice abundance, because strong UV radiation from the early type sources can dissociate 
water molecules located on the mantles of dust grains. This effect is known to arise as the 
threshold extinction {Ay-o) in relationship of N{H20) = q x (Ay - Ay-o) between column 
density of water ice {N[H20)) and visual extinction {Ay). Many studies have been done 
for determination of Ay.Q toward various molecular clouds such as Taurus, p Ophiuchus, R 
CrA, and Serpens (Whittet et al. 1996; Teixeira & Emerson 1999) and Ay.^ is found to vary 
from 2 to 10 for those clouds due to the different internal UV radiation. Teixeira & Emerson 
(1999) found that the q parameter is comparable between two molecular clouds in Taurus 
((1.0 ± 0.1) X 10^^) and p Ophiuhcus ((0.92 ± 0.04) x 10^^), and therefore we adopt g = 1 
X 10^^ hereafter to find the threshold extinction for our targets. 

JW 370 is located at a distance of 246'.'4 from 9^ Ori C with Ay = 5.36, and the water 
ice toward the object is detected with Tjce = 0.12 and Az/ = 443 cm~^ in our data. The 
column density of the water ice is calculated as 2.66 x 10^^ cm~^ using N{H20) = Th^q • 
Az//(2.0 X 10^^^) (d'Hendecourt & Allamandola 1986). As result, the threshold extinction 
Ay.Q for JW 370 is estimated to be 2.70. 

Since dl82-332, dl83-405, and d218-354 are located much closer to the exciting O, B 
stars, much higher threshold extinction should be expected. In this sense, there is little 
chance for dl83-405 and d218-354 to have the water ice mantle around the grains due to the 
lower Ay in these sources (see Table 3). For the case of dl82-332, higher extinction of Ay = 
6.0-6.5 was estimated by McCaughrean & O'Dell (1996). However, dl82-332 is located very 
close to 9^ Ori C (27'.'2) and the much higher threshold of Ay.Q could hamper the water ice 
formation onto the grains. 

d216-0939, which shows the moderate water ice absorption in its disk, has a significantly 
larger distance (400'.'5) from NU Ori (exiting source of M43) compared to the other targets 
in our sample, and it should suffer least from intense UV radiation. Therefore, the threshold 
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extinction of Ay-o can be assumed to be no greater than 2.70. Whereas Hillenbrand (1997) 
derived very small extinction (Ay = 0.72), it is not sufficient to produce the water ice even 
with such low extinction threshold. The Ay of 0.72 is likely to be underestimated due to 
contribution of the reflection nebula in the optical band. 

Lastly, dl21-1925 has moderate distance from 9^ Ori C (247'/l) and the Ay-o should 
be comparable to 2.70 around its vicinity star JW 370. The foreground water ice detected 
toward dl21-1925 is consistent with the estimated Ay of 5.83 (Hillenbrand 1997). Taking 
into account the overall structure of M42 (O'Dell & Harris 2010), the foreground materials 
should act as a good shield from the UV radiation to keep the water ice around the grains 
in the disk. Therefore, the null detection of water ice in the disk of dl2 1-1925 needs to be 
explained in some other way. 

4-2.3. Inclination Angle of Disk 

The most plausible explanation for the water ice features in our targets is related to 
the disk geometry. Pontoppidan et al. (2005), Tcrada ct al. (2007), and Aikawa ct al. (2012) 
clearly detected the water ice features in the disk whose inclination is very high (i > 70°) 
owing to the abundant materials for the absorption. Pontoppidan et al. (2005) pointed out 
that the disk inclination is an important parameter for the water ice detection in their disk 
model of CRBR 2422.8-3423. In their model, the inclination of 69 ± 3° seems to act as the 
critical angle. 

In the disk with an inclination angle smaller than the critical angle, the expected optical 
depth is rapidly decreasing to zero. The two sources: dl21-1925 and dl83-405 are thought 
to have the small inclination angles of 51° and 39°, respectively. No water ice detection in 
the disk of those objects is consistent with the Pontoppidan et al. model. For the remaining 
three sources: dl82-332, d2 18-354, and d2 16-0939 with higher inclination angle of the disks, 
the observed optical depths of the water ice absorption are presented in Figure 4 with an 
overplot of the t^ q^ shown by Pontoppidan et al. (2005). In the case of d2 16-0939 disk 
with the inclination angle of 75-80° above the critical angle, the detected optical depth is 
explained well by this model. 

On the other hand, the question should arise for dl82-332 and d218-354 in this scheme. 
As for dl82-332, the expected optical depth is rather low (< 0.1), and so systematic errors 
for spectral fitting and/or photometric measurements errors could be the cause of this dis- 
crepancy. The inconsistency for the case of d2 18-354 which has rather high inclination angle 
(65°) seems more serious. This should be explained by a higher critical angle of the disk 
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Fig. 4. — Optical depth of the water ice vs. inchnation angle of the disk. Filled square 
with error bar shows the detection of the water ice for the d2 16-0939 disk. Open squares 
represent the estimated upper limit (Icr) for non-detection of the water ice for the dl82-332 
and d218-354. dl21-1925 and dl83-405 with lower inclination angle are not shown in this 
figure. Solid line is the expected optical depth of the water ice in the model for CRBR 
2422-3423 disk derived by Pontoppidan et al. (2005). 
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inclination and/or lower abundance of the water ice. Considering that circumstellar disks so 
far detected exhibit a wide variety of morphology in scattered light, it is possible that the 
same critical angle of the disk inclination as in CRBR 2422.8-3423 should not necessarily be 
applied to all the disks. In fact, the d216-0939 disk seems to have a different critical angle of 
~ 77° judging from its flared disk, because the critical angle is expected to be determined by 
the opening angle of the disk. The problem with both the cases of dl82-332 and d2 18-354 
might be resolved with the optical depth curve shifted accordingly, if they have a higher 
critical angle of the disk inchnation. Regarding the possible difference of the water ice abun- 
dance, we note that it should be affected significantly by much more intense UV field for the 
disks (dl82-332, dl83-405, and d218-354) located at the vicinity of O, B stars. 



4.3. Properties of d216-0939 Disk 

4.3.1. Geometry for Light Path 

As mentioned in the above sections. Smith et al. (2005) revealed that the d216-0939 disk 
has an inclination angle of 75-80° which is well matched to the flared angle of the disk (77°). 
This means that this target might provide with very unique opportunity for us to directly 
measure the absorption feature through the disk toward the central star. In the following, 
this suggested interesting geometry is re-examined using our near infrared data. 

Using the optical depth parameters of Tjce = 0.67 and Au = 538 cm^^, the column 
density of the water ice toward d216-0939 is estimated to be 1.80 x 10^*^ cm~^ in the same 
way as the above. Since the detected water ice is known to reside in the disk region of 
d216-0939, the visual extinction caused by the disk material (Av-disk) along a light- of- sight 
toward d216-0939 can be derived as 20.7 mag. This extinction is usually large enough to 
make the central source undetectable at the optical band. Thus, it is most likely that the 
previously measured optical flux could be considered to originate from the indirect reflected 
light from the central star. In fact, the position of the peak at H and K in our data shows 
a marginal offset into the direction perpendicular to the d2 16-0939 disk, while the peak 
position at K and L' falls into the same area within the measurement accuracy (Figure 5). 
This suggests that the flux shorter than 1.65 ;um appears as the reflected light. However, 
this result apparently conflicts with no detection of the positional shift of the peak flux 
at 0.550 //m and 0.658 /xm bands reported by Smith et al. (2005). Although it might be 
explained by the appropriate size of the large grains, this result is required to be addressed 
by using more sensitive and higher spatial resolution imaging in future. 

Apart from the position offset, the radial profiles of d2 16-0939 at each band are inves- 
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Fig. 5. — Possible change of the peak position of the d2 16-0939 flux. The reference position 
is determined by the position of JW 671: the vicinity star of d216-0939. The peak at the 
longer wavelength shows the possible positional change toward the direction perpendicular 
to the disk. It may suggest that flux at shorter wavelength (at least <1.65 fim) comes from 
the scattered light by material around the star. 
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tigated in Figure 6 to estimate the contribution of the reflected hght to the flux. This flgure 
presents the radial profiles for both directions perpendicular to and along the d216-0939 
disk (position angle: 83° and 173°, respectively). JW 677 (Jones & Walker 1988) which is a 
brighter star located 10'.'84 south is chosen as a reference of the point spread function (PSF). 
In the profiles at H and the extended features are clearly seen for both directions and the 
profiles along the disk appear more extended. Thus, the photometric results at H and K for 
d216-0939 include the scattered light component. On the other hand, at L', no significant 
extended feature can be observed for the direction perpendicular to the disk. The radial 
profile at V along to the disk implies the possible extended flux. Higher signal-to-noise 
imaging with higher spatial resolution is required for concrete conclusion for this extended 
feature. 



4.3.2. Mass & Density 

Considering the lack of change of the K and V peak position (Figure 5) and the lack of 
evident difference in the V proflle of d216-0939 and JW 677 (Figure 6), we assume that the 
L' light reaches us directly from the central star. The derived visual extinction, Av;disk-i and 
column density of the water ice, N{H20), put direct constraints on physical parameters of 
the disk such as mass and density. Firstly, based on the relationship of Nh = (1.59 x 10^^) 
X Ay cm~^ between Nh and Ay (Imanishi et al. 2001), is derived as 3.29 x 10^^ cm~^. 
This number is in good agreement with the relationship between Nh and the inclination 
angle of the silhouette disks in M42 found by Kastner et al. (2005). 

The silhouette disk of d216-0939 exhibits the clear structure of the disk, whose outer 
radius and scale height at the outer edge are measured to be 1'.'3 (572 AU) and 0'.'6 (264 
AU), respectively (Smith et al. 2005). Thus, the path length can be estimated to be 8.56 x 
10^^ cm and the averaged density can be calculated to be 6.38 x 10~^® g cm~^ accordingly. 
If this density is uniformly adopted to the entire disk, the total disk mass is calculated to be 
1.94 x 10~^ Mq. The derived mass is smaller by a factor of 20 than the one from the submm 
measurement. The difference can be naturally explained by the fact that the mid-plane 
of the circumstcllar disk has higher density than the surface of the disk (see for example 
Hogerheijde et al. 2011). 
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Fig. 6. — Radial profiles of d2 16-0939 and JW 677. Upper panels are measured at the 
position angle of 53°to 113°, which is in the direction perpendicular to the d216-0939 disk. 
Lower panels are parallel to the disk direction with the position angle of 143°to 203°. 
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4.4. Nature of Water Ice in d216-0939 Disk 

The most noticeable feature represented in the water ice absorption of the d216-0939 
disk is the peak of the optical depth and its profile. In comparison to the water ice band of 
edge-on disks in Taurus molecular cloud shown in Terada et al. (2007), the overall profile is 
shifted to longward wavelengths together with the peak absorption (Figure 7). In fact, the 
enhancement of the water ice feature around 3.2 /im has never been seen in any astronomical 
sources before, and the cause of this feature could be very pecuhar. The origin of this broad 
absorption in the red part (> 3.4 /xm) could be a scattering effect and/or other moIecTiles such 
as solid methanol and ammonia hydrate (Dartois et al. 2002). However, the determination of 
the red continuum adopted by a black body with a single temperature has some uncertainty 
and it is required to note the possibility that the optical depth at longer wavelengths may 
be overestimated. 

Regarding the peak of the optical depth of the water ice band, it is well known that 
the higher degree of crystallization and larger grain size for water ice shift the optical depth 
peak to a longer wavelength (Smith et al. 1989). The fully crystallized water ice with a grain 
radii of < 0.5 /xm shows a sharper peak around 3.1/im (ex. Dartois & d'Hendecourt 2001). 
There are several sources showing their peak around 3.1 /xm, which is attributed to different 
degrees of crystallization. For the following discussion, it is noted that those sources do 
not show the 3.2 //m absorption. For the case of the d216-0939, the absorption seen around 
3.2 /im is significantly enhanced, whereas the first peak position of the detected water ice 
is located around 3.1 /im which agrees with the signature of the crystallized water ice. The 
peak shift of the optical depth for the fully crystallized water ice was presented by Leger et 
al. (1979). In their work, the fully crystallized water ice with a grain size of 1 /xm shows the 
peak around 3.2 /xm. 

Taking into account fully crystallized water ice with such a large size (~1 /xm), the model 
fitting to the profile of the optical depth of the water ice is considered with the following 
procedure, which is the same method with which Schegerer & Wolf (2010) have obtained a 
sufficient fit to the water ice features found in their sources. The optical depth of the water 
ice absorption is assumed to be a linear composition of the mass extinctions of amorphous, 
and partially and fully crystaUized water ice with a grain size of 0.1, 0.5, 0.8, 1.5 /xm. We 
used the grain model described by Schegerer & Wolf (2010) to calculate the water ice optical 
depth profile, and the mass extinction curves are derived using the program MIEX (Wolf & 
Voshchinnikov 2004) with an assumption of spherical and compact ice particles. Using the 
mass extinction parameters for each type of water ice, the best fit to the optical depth profile 
of the water ice detected in the d2 16-0939 disk is found to result from the combination of the 
large-sized (0.8 /xm) fully crystallized water ice and the small-sized (0.1 /xm) pure amorphous 
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Fig. 7. — Comparison of the normalized optical depths of the water ice between d2 16-0939 
and edge-on disk sources: HK Tau B and HV Tau C. Overall, the optical depth of the water 
ice in the d216-0939 disk is displaced to longer wavelengths. 



- 21 - 



water ice with a mass ratio of 2.5:1. The best fit curve is presented in Figure 8. The grain 
growth indicated by this result is predicted to occur in the protoplanetary disk. Therefore, 
this grain growth signature of the water ice feature strongly supports the origin of the 
detected water ice to be in the disk. As for the grain size of 0.8 fim, it is significantly smaller 
than the grain size (1.9-4 /xm) found by Shuping et al. (2003) for the largest silhouette disk 
dll4-426. It may be due to different environments caused by UV radiation. 

The detection of the fully crystallized water ice indicates that the detected water ice has 
experienced thermal processing to a temperature of ~ 160K. Considering water ice particles 
arc seen with a viewing angle of 75-80° (Smith ct al. 2005), it should be located almost on 
the surface of the flared disk and the temperature of the water ice particles could be heated 
by radiation from the central star. Schegerer & Wolf (2010) interpreted that their spatially 
resolved detection of the crystaUized water ice toward YLW 16A in p Ophuichus molecular 
cloud comes from the outer surface of the disk. In our data, the water ice distribution 
in radial direction cannot be probed directly. Since the grain growth is expected to occur 
most effectively in the mid-plane of the circumstellar disk, it would be natural to consider 
some mixing process to exchange the materials between the mid-plane and the surface of the 
disk. The mid-plane of the circumstellar disk is cold enough to keep water ice in amorphous 
form, and the detected large-size fully crystallized water ice at the surface of the disk likely 
originates from the heating of large-size amorphous water ice transferred from the mid-plane 
of the disk. In this scenario, the detected small-size amorphous water ice particles should 
reside at the colder and less dense outer region of the disk that is further away from the 
central star than the crystallized water ice. 



5. CONCLUSION 

We obtained 1.9-4.2 /xm spectra for five targets associated with the distinct silhouette 
disk. For three objects: dl82-332, dl83-405, and d218-354 no water ice feature was detected 
primarily due to lower angle of the disk inclination than required for detection of the water ice 
signature. In addition those objects are too close to the O, B stars exciting the H II region. 
Whereas dl21-1925 exhibits a slight water ice band at 3.1 fim. with an optical depth of 0.12, 
a nearby star shows almost the same absorption as this object. Thus, the origin of those 
water ice absorption could be common, and it is most likely attributed to the foreground 
materials of the water ice. In contrast, the comparison star near d216-0939 shows no water 
ice feature, while the moderate absorption (r^ce ~ 0.67) of the water ice at 3.1 iim toward the 
d2 16-0939 is definitely seen. Therefore, we concluded that the water ice toward d2 16-0939 
should be localized in ~2000 AU scale and that the material responsible for the absorption 



- 22 - 




Wavelength [ jum ] 

Fig. 8. — Comparison of the normalized optical depth of the water ice between the model 
and the detection. Solid line shows composition of purely crystallized and pure amorphous 
water ice by 2.5:1. This model is found to match sufficiently to explain the absorption profile 
at wavelengths shorter than 3.3 /im. The extra absorption at 3.4 /xm seen in the spectrum 
could be attributed to a scattering effect and/or various possible carriers of the C-H bond 
such as CH3OH ice. 
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could be inside the disk of d216-0939. Due to the clear morphology and suggested ideal 
geometry for our line-of-sight, we directly derived the average density on the disk surface 
to be 6.38 xlO~^^ g cm~^. Remarkably, the profile of the optical depth of the water ice 
shows the signatures of fully crystallized water ice with a grain size of ~ 0.8 fim. Secondary 
composition of small amorphous, that is, non-evolved water ice indicates that this material 
could exist in outer region of the circumstellar disk. 

The water ice feature at 3.1/xm is at a wavelength where diffraction-limited images 
can be easily obtained by using adaptive optics (AO) with the 8-m Subaru telescope. The 
spatial resolution of ~ 0'.'08 at 3.1 fim is comparable to the resolution achieved by HST in 
the optical. Future AO observations of the silhouette disks can help to reveal the morphology 
of the water ice absorption in more detail. 

We thank all the Subaru telescope staff for maintaining the IRCS and supporting the 
observation. We also acknowledge the development efforts of the instrument to the develop- 
ment team members to realize these high sensitivity observations. 
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Table 2. Magnitudes of Targets 



Object 


H 


K 


L' 




(mag) 


(mag) 


(mag) 


dl21-1925 




11.74^ 


10.83 


dl82-316 


13.10^ 


12.25^ 


11.18 


dl83-405 


12.80^ 


12.24^ 


10.93 


d218-354 


9.69^ 


9.49^ 


9.00 


d216-0939 


13.29 


12.49 


10.88 



^Muench et al. (2002). 
•^McCaughrean & O'Dell (1996). 



Note. — The typical 1 a uncertain- 
ties are 0.02, 0.02, and 0.1 mag for 
i?, and L' in our data, respectively. 
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Table 3. Spectral Type, Visual Extinction, Blackbody Temperature and Infrared Excess 





Previous Works 




Derived Parameter 






Spectral 




Spectral 




Blackbody 




Object 


Type 


Av 


Type 


Av 


Temperature 


rL' 






(mag) 




(mag) 


(K) 




dl21-1925 


M4.5^ 




M6 


6.8 







dl82-332 




e.o-e.s'^ 


M6 


6.5 


1000 


0.13 


dl83-405 


M3^ 


2.55^ 


M3 


2.55 


1100 


0.29 


d216-0939 




0.72'' 


K5 


0.72 


1100 


1.7 


d218-354 


G6-K3'^ 


1.51^ 


Ml 


1.3 


500 


0.29 


Vicinity Stars 


JW 370 


K0-K3'' 


5.36'' 


Kl 


5.36 


1000 


0.19 


JW 071 


AI2'' 


1.72'' 


M6 


1.72 


l.")00 


0.14 



'^Getman et al. (2005). 
''Hillenbrand (1997). 
^McCaughrean & O'Dell (1996). 

Note. — r^/ is defined as the ratio of the infrared excess at 3.77/im to the 
intrinsic stellar continuum flux: Fe2;cess(3.77/zm)/Fj„tj.i„ric(3.77/xm). 



Table 4. Optical Depth of Water Ice Absorption at 3.1 //m 



Object 


Peak 




Az/ 
(cm~^) 


dl21-1925 


3.10 


0.12±0.01 


443 


dl82-316 




<0.012 




dl83-405 




<0.022 




d2 18-354 




<0.054 




d2 16-0939 


3.12 


0.67±0.01 


538 
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Table 5. Disk Parameters 









Inclination 


Exciting 




Object 


Mass*^ 


Diameter 


Angle 


Source 


Distance 




(xlO-2 Mq) 


(") 


n 




(") 


dl21-1925 


0.74±0.07 




51'' 


6*1 Ori C 


247.08 


dl82-332 


<0.37 


0.3^ 


55'' 


Ori C 


27.23 


dl83-405 


<0.20 


0.7'' 


39'' 


Ori C 


50.28 


d216-0939 


4.53±0.06 


2.6'= 


75-80'= 


NU Ori 


400.49 


d218-354 


2.37±0.04 


1.4'' 


65'' 


6*1 Ori C 


85.39 



"^Mann & Williams (2010). 
''Bally et al. (2000). 
'=Smith et al. (2005). 



Note. — Distance is measured from the listed O, B stars responsible for 
excitation of the H II region. 



